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We prove the existence of an off-diagonal long-range order in a one- 
■ dimensional many-body problem. 
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It is well known that in a many body system of bosons or fermions it 
is possible to have an off-diagonal long-range order (ODLRO) of the re- 
duced density matrices in coordinate space representation Q. The onset 
of ODLRO points at quantum phases and quantum phase transitions in a 
many-body system. For a many-body bosonic system, the well-known ex- 
ample is the Bose-Einstein condensation (BEC). It is reasonable to assume 
that superfluid He II and superconductors are quantum phases character- 
ized by the existence of such an order. A generalized criterion developed by 
Penrose and Onsager states that a system exhibits ODLRO if the largest 
eigenvalue of the one-particle reduced density matrix is extensive Q . Besides 
these well-known examples of three-dimensional systems, it was relatively 
recently shown that there exists a novel type of ODLRO in the fractional 
quantum Hall effect ground state ||]. Related to this, presence of ODLRO 
in the ground state of a Calogero-type model in two dimensions has also 
been demonstrated @. 

Recently, two of us (JK) || have presented a many-body system in one 
dimension with upto next-to-nearest neighbour interaction and obtained its 
exact ground state. For the version of the model where N particles reside 
on a circle, they calculated the one and two point correlation functions 
by using this ground state and showed that there is no long range order. 
Further, by applying the Penrose-Onsager criterion it was claimed that there 
is no ODLRO as well. Unfortunately, the Hamiltonian employed there to 
describe the system is not a symmetrised one and since the Penrose-Onsager 
criterion does not apply for the case of distinguishable particles, the claim 
regarding the absence of ODLRO is incorrect. In this Letter, we present the 
symmetrised version of the many-body system and prove that indeed there 
exists ODLRO. To the best of our knowledge, this is the first example of a 
one-dimensional system to possess ODLRO and hence, quantum phases. 

Consider N particles on a circle, described by the symmetrised Hamil- 
tonian, 

1 N d 2 
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where O is the step function and W{x\, ...,x^) is the iV-body potential 
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In (|), the sum is over all permutations on N symbols in the symmetric 
group, Sn- The exact ground state of the un-symmetrised version of (|l|) 
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was obtained in ||. Relying on the ground state found there, we introduce 
the symmetrized (un-normalized) wave function : 

1pN(xi, ...,X N ) = (p N (x P{ i ) , ...,Xp( N )) , (3) 

where P is the permutation in Sjy such that 1 > Xp(i) > ip( 2 ) > ... > 
X P(N) > 0) an d 4>n is the (un-symmetrized) exact ground state wave func- 
tion of the un-normalized Hamiltonian which is given by 

TV 

(f> N (xi,...,x N ) = Y[ I sin7r(x n - x n +t) f ] (x N+1 = x±) , (4) 

n=l 

provided g = (3{f3 — 1), G = (3 2 (we have set the scale factor L equal to 1). 
Primitively, the function (||) is defined on the hypercube [0, 1]^. 

Now for (3 > 2 it is easily verified that -0/v can be continued to a multi- 
periodic function in the whole space 7Z N (or equivalently on the torus T ): 

ip N (x 1 ,...,x i + l,...,x N ) = ip N (xx,...,Xi,...,x N ); (i = 1, N) , (5) 

which belongs to C 2 (i.e. is twice continuously differentiable) . Owing to this 
property and the results of ||, ipw then obeys the Schrodinger equation (with 
Hamiltonian (1) and with the same ground state energy Eq = N((3tt) 2 ) not 
only in the sector x\ > X2 > ... > xjy but everywhere. Thus, tp^ describes 
the ground state wave function of the iV-boson system. Moreover, it is 
translation- invariant (on 1Z N ): 

iPn{xi + a,x 2 + a, ...,xn + a) = iPn(xi,x 2 , x N ) ; V a G TZ . (6) 

We are interested in the one-particle reduced density matrix, given by 

p N (x-x') = — I dx\... \ dx N -i 

JO JO 

-0/v(xi, ...,XN-l,x)lf) N (xi, ...,x N _i,x') , (7) 
where Cjy stands for the squared norm of the wave function: 

Cjv = / dxi... [ dx N \ i(j n (xi,...,x n ) \ 2 . (8) 
Jo Jo 

That the right hand side (RHS) of Eq. (0) defines a (periodic) function of 
(x — x') is an easy consequence of Eqs. (||) and (|6|). The normalization of pn 
is such that /9jv(0) = N, the particle density. Further, the function pn(0 is 
manifestly of positive type on the U(l) group, which implies that its Fourier 
coefficients, 

pP = f 1 d£e- 2i ™Sp N (0 ; (n = 0, ±1, ±2, ...) , (9) 
Jo 
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are non-negative (Bochner's theorem). In fact, this directly appears if one 
writes their explicit expression in the form (obtained by using the periodicity 
property): 



(n)_ N 

Pn — r< 
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(10) 



Since the function p^ is not only of positive type but also positive (like 
-i/V), Eq. (|9[) shows us that 

pP>pPl (n = ±l,±2,...). (11) 

(n) 

Notice that the coefficients p N , which physically represent the expec- 
tation values of the number of particles having momentum k n = 2im in 
the ground state, are nothing but the eigenvalues of the one-particle re- 
duced density matrix (diagonal in the k n representation). According to the 
Penrose-Onsager criterion 0], no condensation can occur in the system (at 
least for Bose particles) if the largest of these eigenvalues is not an extensive 
quantity in the thermodynamic limit, that is, if 

(0) 

l im ^L = o. (12) 
Accordingly, we have to evaluate 



N C N ' 
where CV is given by Eq. (^) and 



(13) 



A N = [ dxi... [ dx N -iip N (xi, ...,x N -i,0) 
Jo Jo 

fl 

dxip N (xi, ...,x N -!,x). (14) 



Because of the special form, of the wave function, the computa- 

tion of the squared norm Cat is more difficult than in the case of iV free, 
impenetrable particles. As a consequence, the (mainly algebraic) method 
introduced long ago by Lenard [26] to deal with the latter case does not 
apply here, and we have to resort to another device. For conciseness, we 
introduce the notation: 

S(x n - x n+ i)= | sin7r(x„ - x n+l ) f , (15) 

and define: 

,A 

5 2 (A) = / dxS(x)S{A -x); (0 < A < 1) . (16) 
Jo 
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Our starting point will be the following representations of Cm and A jy: 



C N = {N-l)\ — 



dxe lx F{x 
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(N-iy. 
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dxe lx F(x 



,N-3 



F(x)G{x)+(N-2)H(xY 



(17) 



,(18) 



where 



F(x) 
G(x) 
H(x) 



dAe iAx S(A) 



dAe iAx S 2 (Af , 
dAe iAx S(A)S 2 (A) . 



(19) 



The representations (|l7| )-(|i"9|) follow from the convolution structure of the 
expressions (§) and ([lj) of Cm and Am, when written in terms of appropriate 
variables. Their proof will be given elsewhere Q. Our aim is to extract from 
them the large N behaviour of Cm and A m- Their form is especially suited 
for this purpose, because the integrands in Eqs. (|l7| ) and (|lg|) are entire 
functions, as polynomial combinations of Fourier transforms of functions 
with compact support (Eq. (|i~9|)). Indeed, we are then allowed to, first, 
shift the integrand and then apply the residue theorem to meromorphic 
pieces of the integrands. However, it turns out that the calculations needed 
for arbitrary (integer) values of (5 are quite cumbersome. So, in order to 
keep the argument clear enough, we shall content ourselves to present below 
these calculations in the simplest case, namely (3 = 1 (recall that, strictly 
speaking, this value is not allowed), being understood that similar results 
are obtained for all integers (3 > 2. To be sure, after the illustration of 
the calculation for the case of j3 = 1, we give the final results for arbitrary 
integer value of (5. 

For (3 = 1, S(A) = sin7rA, and Eq. (|l9| ) gives, after reductions: 



F(x) 
G(x) 
H{x) 



2tt 2 



l x\x c 
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4vr 3 



i x(x 2 — 47T 2 ) 2 



+ e' 



^(" 2 )| 



X 



(20) 



where R^ n \x) is a generic notation for rational functions behaving like x n 
when x — > oo, and the precise form of which will be eventually of no impor- 
tance. This produces, for the functions to be integrated in Eqs. ([T^ ) and 
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2tt 



2' 



N 



F(x 



x(x 2 -Att 2 )} n _ l 
F(x) N - 3 [F(x)G(x) + (N- 2)H{xf] 

^2\ N rfOAT , 1U AJ* N+l 



1 N 

+ V e ina: i2 ( - 37V) fa;) 



(21) 



Let us stress again that these functions, when analytically continued, are 
holomorphic in the whole complex plane (the poles appearing in the first 
term are exactly cancelled by the remaining ones). 
We consider first Cat, now given by 



c N = (N-iy. 



2tt 2 \ 



oo N 
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27T J-oo I [x(x 



dxe ~ lX{ ' ' " \ 7T 2 )] N + E e^T^O*)} • (23) 



Since the function within the curly bracket is an entire one, we can shift the 
integration path to I = {z = x + ia \ x e TZ}. Let us choose a > 0. Then, 
by Cauchy theorem 

r N 

/ dze~ lz e inz Ri~ 3N) (z)=0. (24) 

^ n=l 

Indeed, the integrand is holomorphic above / and is bounded there by const. 
| z \~ 3N , which allows us to close the integration path at infinity in the 
upper complex plane. We end up with 

/2tt 2 \ N 1 r e~ iz 
C K=(N -!),(_] _^ 2 _ w _ 7 _ w . (25 ) 

Similarly, we are allowed to close the integration path at infinity in Eq. 
(25), but this time in the lower complex plane. For obtaining the large- N 
asymptotics, we write 

where R±(a) and Ro(a) are the residues of the last integrand at z = ±y/a 
and z = respectively. They are readily computed and sum up to 

R+(a) + + R (a) = (-1) M+1 V - -a s (27) 

+ \ J \ ) uv ) v ) (2M + 2s + 2 ! v ; 
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for N = 2M + 1. 

Using Eqs. (H), @ and @ we then obtain 

JV 




7 a 



n t^ n!(3iV + 2n-l)! 



(2M + 2s + 2)! 
2 ^ n (28) 



The result is exactly the same for even N. It suffices now to observe that 
the last series alternates in sign and is decreasing to deduce 



C N = (2tH) 



2WV (^"1)1 



(3iV-l)! 



l + 0( 



u 
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(29) 



Our procedure for evaluating An is quite similar, and we give below the 
final result, viz., 



A N = (2N+ l)(2vr 



2^(^-1)1 



(3N)\ 



1+ < 



Finally, using Eqs. (O), (|29|) and (pO), we obtain 



(0) 

Pn 
N 



1 + 



(30) 



(31) 



The same procedure applies for all integer values of f3, although the algebra 
becomes quite involved. The general result (for any integer (5) is: 



lim 

TV— >oo 



(0) 

Pn 



(/3!) 4 [(3/3 + l)!] 2 



(32) 



N [(2/3)!]2[(2/? + l)!]3- 

Our method does not generalize in a straightforward manner to the case of 
non-integer values of /3, but there is clearly no reason to expect a different 
outcome for such intermediate values. 

To summarize, the Penrose-Onsager criterion di~2] ) is not met for bosons 
in this model as the limit of the ratio in (^) is non-zero. The onset of 
such an order leads to a new thermodynamic phase of the system. Thus, we 
reach the remarkable conclusion that Bose-Einstein condensation is indeed 
possible in the bosonic version of the iV-body model discussed in ||. We 
recall once again that, to the best of our knowledge, this is the only example 
in one-dimensional statistical mechanics where there exists ODLRO and 
hence quantum phases. 
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